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66 ■ 

(^— N , We study the properties of the ground state of Nonlinear Schrodinger Equations with spatially 

^— V ■ inhomogeneous interactions and show that it experiences a strong localization on the spatial region 

fvj ' where the interactions vanish. At the same time, tunneling to regions with positive values of the 

interactions is strongly supressed by the nonlinear interactions and as the number of particles is 

^ . increased it saturates in the region of finite interaction values. The chemical potential has a cutoff 

■ ' value in these systems and thus takes values on a finite interval. The applicability of the phenomenon 

to Bose-Einstein condensates is discussed in detail. 
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c/^ '. i. introduction 

Oh 

^ ■ The experimental generation of Bose-Einstein condensates (BEG) with ultracold dilute atomic vapors [l| has turned 

out to be of exceptional importance for physics. The formation of a condensate occurs when the temperature is low 

^ C^ ■ enough and most of the atoms occupy the ground state of the system. This process is visible both in momentum and 

I in real space due to the spatial inhomogeneities exhibited by the order parameter on a macroscopic scale because of 

y—{ , the trapping potentials. 

^ ' The properties of the ground state of trapped BECs are well known. In the mean field limit simple analytical 

^^ , expressions are available in the Thomas-Fermi approximation 2| and beyond [3, 4] . These approximations and direct 
^^ • numerical simulations describe accurately the properties of the experimentally found ground states [5|, |6| . 
^-j ' Nonlinear interactions between atoms in a Bose-Einstein condensate are dominated by the two-body collisions that 

. ' can be controlled by the so-called Feschbach resonance (FR) management Ml. The control in time of the scattering 
^T ' length has been used to generate bright solitons @, |^, Upland induce collapse fill and has been the basis for theoretical 
^^ . proposals to obtain different types of nonlinear waves [i3,[ij, H, [laMj[l3j[lM- Moreover, interactions can be made 
spatially dependent by acting on either the spatial dependence of the magnetic field or the laser intensity (in the case 
of optical control of FRs [l^) which act on the Feschbach resonances. This possibility has motivated many theoretical 
studies on the behavior of solitons in Bose-Einstein condensates (BECs) with spatially inhomogeneous interactions 

[11 [IE [13, 111, [11, 111, [11, [13, IM 111, H ill, 111, [11, 111 . 

5J] ■ In this paper we describe a novel nonlinear phenomenon related to the ground state of a BEC in the mean field 

Cd ' limit when the interactions are spatially inhomogeneous. We will show that when the scattering length is non-negative 
and vanishes on a certain spatial region, a striking localization phenomenon of the atom density occurs at the regions 
where the interactions vanish. By tuning appropriately the control (magnetic or optical) fields this phenomenon can 
be used to design regions with large particle densities and prescribed geometries. Another interesting phenomenon to 
be studied in this letter is the nonlinear limitation of tunneling of atoms to the regions in which the interactions are 
stronger and the fact that the chemical potential exhibits a cutoff value. 

Although it is reasonable to think that a BEC will avoid regions of large repulsive interactions and prefer to remain 
in regions with smaller interactions, the localization phenomenon to be described in this paper here goes beyond what 
one would naively expect. 



II. THE PROBLEM AND ITS MATHEMATICAL MODELLING 
A. Physical system to be studied 

In this paper we will consider physical systems ruled by the nonlinear Schrodinger equation (NLS) 

i^ = -lA^ + V{x)^P + gixMP^, (1) 

in M.^ , where p > 1 is a real parameter and g, V > arc non-negative real functions. V describes an external localized 
potential acting on the system satisfying, 

V{x) -^ oo, as |x| -^ oo, (2) 

and 5 is a spatially dependent nonlinear coefficient. Stationary solutions of Eq. ([1]) are defined through 

ip{x, t) — (j){x) exp (lAt) (3) 

which leads to 

X<f>^^^A^ + V{x)c^ + g{xM\P^. (4) 

Of all the posible solutions of Eq. ([4]) we will be interested on the so-called ground state, which is the the real, 
stationary positive solution of the Gross-Pitaevskii equation (j4]) which minimizes the energy functional 



Ei4>) 



l\V^\^ + V{x)\ct,\' + ^gix)\<t>r^ 



(5) 



under the constraint of a fixed J^3 \(j)\'^dx. 

For our purpouses the only relevant property is the positivity, thus we will consider positive solutions of 

-^Au + V{x)u^Xu-g{x)uP, x e R^ (6) 

where V satisfies ([2]). 

From the physical point of view Eq. ([5]) is a model of different physical systems. However, as stated in the 
introduction, in this paper we will focus mainly on its applicability to the description of the ground state of a trapped 
Bose-Einstein condensate with spatially inhomogeneous interactions in the mean field limit where physically it is both 
posible to control the interactions and achieve large values for the coefficient g{x). Although Eq. (HI) is written in 
nondimensional units it is easy to connect this problem with the realistic quantum problem (see e.g. [l|), considering 
that the coordinates x and time t are measured in units of ao = \/h/rrvw and 1/w, respectively, while the energies 
and frequencies are measured in units of hw and w respectively, w being a characteristic frequency of the potential. 
Finally, g{x) = 47ras(x)/ao is proportional to the local value of the s-wave scattering length as{x) and the parameter 
p = 2. Another relevant physical quantity, the number of particles, is directly related to the L^ norm of the solution 
through 

N=\\u\\l^ f u^dx. (7) 

In this paper we will study properties of the positive solution of Eq. ^ (ground state) when the interactions vanish 
on a certain set of points which we will denote as w, i.e. 

w := {x e R^ : g{x) ^ 0}, (8) 

Through this paper we will assume that lo is composed by a finite number of closed connected components ldj^ 
1 < j < J, which are disjoint LUi n loj = if i 7^ j and it is assumed that each component loj is regular enough and 
that R^ \ oj is connected for iV > 1. 



The Thomas-Fermi limit 



The Thomas-Fermi approximation is used extensively in the field of Bose-Einstein condensation to describe the 
solutions in the situation in which the nonlinear term is "large" and proceeds by neglecting the kinetic energy or 
equivalently the term proportional to A(/) in Eq. ^ which leads to 



4>TFix) 



l\~V{x) 
9{x) 



(9) 



Thus, the TF approximation provides the profile of the stationary positive solution (i.e. the ground state). In the 
case of spatially homogeneous interactions g{x) = ga the so-called classical turning points for which A = V{x) delimit 
the regions in which the approximation breaks, since near to those points the amplitude is small and thus the kinetic 
energy cannot be neglected anymore. 

It is remarkable that, when the interactions are spatially inhomogeneous a new feature appears, which is that 
the amplitude of the Thomas-Fermi solution diverges on the set w. Although it is obvious, that in the vicinity of 
Lu the Thomas- Fermi approximation is not correct, this divergence is the first indication of an interesting nonlinear 
phenomenon with relevant physical implications whose study is the purpouse of this paper: the tendency of positive 
solutions of Eq. (0j, or in physical terms the ground state of Bose-Einstein condensates with spatially inhomogeneous 
interactions, to localize strongly on the regions where the scattering length is close to zero provided the system is 
sufficiently nonlinear. 



III. LOCALIZATION PHENOMENA: A "TOY" EXAMPLE 



Let us first consider the exactly solvable one dimensional "toy" model on the interval x G [—L, L], given by 

Xu = —u" + g{x)u , 
u{L) = u(-L) = 0, 

which can be understood as a version of Eq. © with a potential of the form 

V{x)^ 
We will take the interactions as given by the equation 

9{x) = 




go \x\ < a, 
\x\ > a. 



(10a) 
(10b) 



(11) 



(12) 



This is a model for an ideal quasi-one dimensional BEC in a box. 

In this simple case, the positive, stationary solution, can be obtained analitically and is given by 



,a < \x\ < L, 



u{x) — C sin V2X{x ~ L) 
while for |a;| < a the solutions are given by 



(13a) 



^ sn(xVXa + S;k' 
uix)^{7r/{2,/2^\a-L\), 

I xvXa; k^ j , 



'^ dc ( a;\/Aa; fc2 ' 



A < A, 

A ^ A:^; 

A > A* 



(13b) 



where a{k) = 2/(1 -I- fc^), sn and dc are two of the standard Jacobi elliptic functions and k is the elliptic modulus. 
Both the elliptic modulus and amplitude C can be obtained from the matching conditions for ^(a) and (f>'{a). These 
conditions also give A, = tt^/ [8(a — L)^] . 

It is interesting to point out that there is a cutoff value for the eigenvalue, Ac beyond which stationary solutions 
do not exist. Its explicit value can be obtained from the condition of maximum slope at a; = a which leads to 




FIG. 1: [Color online] Spatial distribution of n{x) = u (x) for the ground state solutions given by Eqs. Eq. (I13ap - (|13b[) with 
L = 2,a — 1 and different values of A corresponding to different values of A^ (a) From the lower to the upper curve goN = 2, 3.55 
(corresponding to A.) and goN = 6. (b) From the lower to the upper curve goN = 25, 50, 100, 200. 



Ac = 7r^/[2(a — L)'^]. In that limit the i^-nomi (number of particles in the condensate) diverges, since the amplitude 
in the outer region C — > oo. 

The spatial profiles of the ground state density for different values of goN shown in Fig. [T] support our conjecture 
based on the Thomas-Fermi solution, i.e. the existence of a strong localization of the solution in the region where the 
interactions vanish. 

Physically it is also remarkable that the atom density in the inner part of the domain, i.e. the region where there 




FIG. 2: [Color online] (a) Ratio between the maximum atom density and the atom density at a; = as a function of the scaled 
scattering length g^ on the spatial region |a;| > a (for go = 1) (b-c) Scaled atom density profiles for (b) g* = and (c) g, = 0.2 
in both cases for a total scaled number of particles A'' = 1000. The insets show the profile of g(x). 

are nonlinear interactions, remains almost constant independently on the number of particles in the condensate once 
a certain critical density is achieved. This region is a nonlinear analogue of the classically forbidden region in ordinary 
potentials and is energetically less favourable due to the extra repulsive energy provided by the nonlinear interactions. 
However, the tunneling of atoms in this region is essentially limited to a constant value, independently of the number 
of atoms, which differs essentially from ordinary tunneling. 

The supression of tunneling strongly depends on the value of the scattering length in the outer region, that we 
will take to be nonzero in what follows. If instead we set g{x) = 5, when |a;| > a and study the dependence of the 
ratio between the maximum atom density and the atom density at a; = (which is a measure of the amplitude of the 
tunneling), we find a strong dependence on this parameter as it is shown in Fig. [Ua). 

This effect is also seen in the atom density profiles when comparing the cases with g* —Q [Fig. [21(b)] and g* = 0.2 
[Fig. EJc)] for N = 1000. Larger values of N lead to a stronger effect. 

Thus we have shown a simple example which already presents the main features of the nonlinear phenomena to 
be studied with more detail in this paper: the localization of the ground state in the regions where the interactions 
vanish, the supression of tunneling induced by the presence of inhomogeneous interactions and the limited range of 
variation for the eigenvalue A (chemical potential). 

Related phenomena have been described in the mathematical analysis of similar logistic equations [35| ] in bounded 
domains [36] . However to consider situations of real physical interest we must move to unbounded domains where the 
analysis is much more complicated and goes beyond our particular simple example and previous theoretical studies 



IV. NUMERICAL RESULTS 



A. Introduction 



In this section we will consider more realistic situations corresponding to problems posed on unbounded domains 
and strongly confining potentials satisfying the condition ^. In the examples to be presented in this section we 
will use harmonic potentials V{x) ex x^ which arise naturally in the modelling of magnetic trapping potentials for 
ultracold atoms. 

To study the localization phenomenon we have computed numerically the ground state solutions of equations of 
the form ([1]) using a standard steepest descent algorithm on the energy functional complemented with a projection 
over the set of solutions of given L^-norm [37J . 
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FIG. 3: [Color online] Ground states of Eq. (g]) for V{x) = 0.02x^ and go{x) = l,gi{x) = exp(-2:^/200), and g2{x) = 
exp(— x^/50), for different values of the scaled number of particles N. (a) Maximum particle density maxa, n{x) and (b) width 
W = J X n{x)dx/N for go (green), gi{x) (red) and g2{x) (blue), (c) Spatial profiles of ^n{x) for A'^ = 10 (blue), N — 100 
(green). A'' = 1000 (red) for g(x) = gi(x). The dashed black line is the ground state with homogeneous interactions and 
A^ = 1000. (d) Same as (c) but for N = 10000 (blue), and A'^ = 40000 (green), in comparison with the case of spatially 
homogeneous interactions (dashed black lines). 

We will consider nonlinear coefRcients g{x) also of realistic forms which vanish or are very small on certain sets. In 
this section we will use a more physical language thinking on the applications of our results to BEC systems. 



B. One dimensional systems 



Let us consider Eq. ([4]) in one-dimensional scenarios [N — 1), and take the potential to be of the form V{x) — 0.02a;^. 
Because of the possibilities for optical control of nonlinearities which are posible in BEC systems we will choose the 
nonlinear coefficient to be of the form 



ga{x) = exp(-a;^/2a^) 



(14) 



for different values of a < -l-oo. These choices will allow us to study different degrees of localization of the interactions 
starting from the case of no localization. 

Our results are summarized in Fig. [31 First, in Fig. [3I^a) we observe how the maximum density (n{x) = m^(x)) 
increases drastically for spatially decaying nonlinearities (blue and red curves) as a function of the effective number 
of particles in the quasi-one dimensional condensate N . This amplitude growth is due to a strong localization effect 
near the region where g{x) vanishes as it is seen in Fig. [31[c,d). In contrast, the condensate density for spatially 
homogeneous interactions grows slowly according to the Thomas- Fermi prediction maxn(x) ex iV^". When the 
number of particles is small, the size of the atomic cloud is smaller than the localization region of g{x). For larger 
values of N the ground state extends beyond the localization region of g{x) and the atom density becomes more 
and more localized near its edge. This effect is more clear for larger N and is accompanied by a saturation in the 
amplitude growth in the region where g{x) is far from zero [Fig. E^d)]. 

In the case of spatially homogeneous interactions the width grows according to the law W ex 7V^/^ [Fig. Efb)]. 
When interactions are spatially dependent and due to the localization of the amplitude close to the zero of g{x) the 
width growth saturates for large values of A^ to a value depending on the size of g{x). 

It is interesting to point out that although the nonlinearity does not vanish strictly anywhere there is an effect 
similar to the existence of the cutoff discussed in the previous example, the number of particles increasing enormously 
for increasing values of A, as it is shown in Fig. 3] 



200000 



160000 



TV 



120000 



80000 



40000 



1 r\ 




/ (a) 


1.0 
0.5 

- °. 


(b) / 


) 0.1 0.2 

A 





10 20 


30 


40 



A 



FIG. 4: [Color online] (a-b) Dependence of the number of particles in the ground state A'^ on the eigenvalue (chemical potential) 



A for g{x) = exp(- 
the linear limit. 



-X /200), and V{x) = 0.02a; . Subplot (b) shows the behavior in for small A'' where the eigenvalue approaches 



These effects are even more clear when the nonlincarity decays to zero faster and is exactly zero on a certain region. 
For instance, taking a nonlinear coefficient given by (74(2;) = (1 — 0.001x^)_|_ (i.e. an inverted parabola with maximum 
amplitude g = 1 at x = and zero values for |a;| > 31.6) as shown in Fig. [5l we see that the maximum densities 
are even higher than for gi(x), g2{x). In this case, in comparison with our first simple example, we can see an even 
stronger localization since the existence of the potential makes energetically more favourable the localization close to 
the point where the interactions vanish. 

In this case we also observe the existence of a cutoff in the values of the eigenvalue A (chemical potential) close to 
a certain value Ac. An example is shown in Fig. [6l 

As a final example, in Fig. [7] we show again the localization phenomenon but now near a; = for localized 
interactions given by g{x) = 1 — exp(— x^/SO) and V{x) = 0.02a;^. In this case interactions vanish in a single point 
a; = being very small in its vicinity. It is interesting to see how localized the density becomes to "avoid" penetrating 
into the regions with appreciable values of the scattering length. 



C. A three-dimensional example 



The same localization phenomenon is also present in three-dimensional geometries. 

As a simple example we will take a set of parameters of applicability to BECs describing a ^^Rb condensate where 
the nonlinear interactions gN are typically in the range 10^ — 10^ for typical numbers of atoms and trap sizes. As to 
the potential and nonlincarity we will take them to be of the form 



1 



V{x,y,z) = -(x^ + j/ + z^) 



g{x,y,z) = go 



1 — exp 



x"^ + y"^ + z^ 



2u;2 



(15a) 
(15b) 



with w ~ 0.5 (i.e. half the trap width). Our results are summarized in Fig. [5] where again the same localization 
phenomenon, here inside the "hole" of the nonlinear coefficient, is seen. For larger number of atoms the localization 
phenomenon can be even enhanced as in the previous one dimensional examples. 
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FIG. 5: [Color online] Ground state maximum particle density (a) and width W = f x^n{x)dx/N (b) for V{x) = 0.02a::'^ and 
go{x) = 1 (red lines) and g4{x) = (1 — 0.001a;^)+ (blue lines). 



V. THEORY 



Introduction 



In Sections IIIII and IIVI we have presented an analytic simple example and several numerical simulations showing 
a striking localization of the solutions of Eq. ^ close to the zeros of the nonlinear coefficient g{x). Here, we will 
provide a rigorous support for two of the observed phenomena: (i) the existence of a finite range of values of the 
chemical potential (i.e. the existence of a cutoff) and (ii) the divergence of the amplitude of the solution close to that 
point. 

In this section we will use a more formal and precise language and provide some mathematically rigorous results 
which will support our results in general scenarios. 
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FIG. 6: [Color online] Dependence of the number of particles in the ground state A^ on the eigenvalue (chemical potential) A 
for g4{x) = (1 - 0.0011^)+, and Vix) = 0.02a;^ 



B. Preliminaries and notation 



We shall fix the potential V satisfying hypothesis ^. Following [3g|, let Q C M^ be an open nonempty set, possibly 
unbounded, with boundary regular enough, let us denote by 



H{n,V):^ lueH^in) : f V u"^ < +oo\ . 



H{il., V) is the completion of C^{n) in the metric derived from the norm 

1/2 



\u\\ := I / \\/u\^ + Vu' 
In 



and H{n, V) is a Hilbert space with the scalar product 



(16) 

(17) 

(18) 

(19) 
(20) 



{u,v) := / \7u\/v + Vuv. 

Let us consider the operator L defined by 

Lu := -Au + Vu, ior ue D{L,n) 
where 

D{L,n):=H{n,v)n{ueL^{n) -. -Au + Vu ei'^in)}. 

Let us first write the following lemma 

Lemma 1. Let 51 C M.^ be an open nonempty set, possibly unbounded, with boundary regular enough. If V satisfy 
hypothesis ^, then the following assertions are true 

i) H{fl,V) ^^ L'^{fl) with compact embedding 

ii) the operator L defined by Eq. il9\) . has a discrete spectrum, noted by cr(L,r2), i.e. a(L,n) consists of an infinite 
sequence of isolated eigenvalues {an{^)} T oo with finite multiplicities. 
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FIG. 7; [Color online] Ground states for V{x) — 0.02a::'^ and g{x) = 1 — exp(— a;^/50) (shown in panel (c)). (a) ^n{x) for 
N = 10 (blue line), A'^ — 100 (green line) and A'^ = 1000 (red line), (b) Condensate width as a function of A'^. 



iii) Moreover, the Rayleigh sup-inf characterization for the eigenvalues holds, and in particular the first eigenvalue, 
denoted by ai{fl), satisfies 



ai{n) 



inf 

tp£H(n.y) 



iVi/jl^ + Vip^ 



4>' 



(21) 



iv) the first eigenvalue is positive, simple with a positive eigenf unction, 0i(f2) > 0, and there is no other eigenvalue 
with a positive eigenfunction. 

This lemma consists of known results, being folk for A^ = 1,2. In higher dimensions, parts i)-iii) are contained in 
[38l | for N > 2, more precisely, see [38|, Theorem 1.2] and the beginning of its proof for parts i) and ii) respectively, 
see [Sa, Lemma 2.1] for part iii). With respect to part iv), it is also well known that the elliptic operator L as defined 
in (|19p admits a unique principal eigenvalue in Q. This principal eigenvalue is the bottom of the spectrum of L in 
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FIG. 8: [Color online] (a) Comparison of the profile of yn(aJ7o7o) for goN — 10^ (blue line) and goN — 10* (red line) 
showing the localization of the atom density, (b-c) Pseudocolor plot of y/n{x,y, 0) for (b) goN = 10^ on the spatial region 
ix,y) e [-10, 10] X [-10, 10] and (c) goN = W^ on the spatial region {x,y) G [-15, 15] x [-15, 15] 



the functional space H{Q, V), and it admits an associated positive principal eigenfunction, denoted by 0i > 0. This 
result follows from the Krein-Rutman theory and from compactness arguments (cf. [331 ^^'^ EOl)- 

In general, the domain of the operator, the eigenvalues and the eigenfunctions depend not only on the domain Q but 
also on the potential V i.e. D{L, H.) = D{L, H., V), cr„(ri) — CT„(ri, V), (/)„(r2) = (j)n{^, V) and so on. We will consider 
a fixed V and in what follows we will skip the dependence on V to clarify the notation, except for the ambiguous 
cases. 

The following lemma is the Lax-Milgram lemma, see [4l|, Lemma V.8]. Let [H{n, V)]* denote the dual space of all 
linear and continuous functionals defined on -ff (il, V), 

Lemma 2. For any f G [H{Q, V)] there exists a unique u G H[U,, V) such that 



n Jn 



(22) 



Remark 1. The above Lemma can be understood in the following way, the inverse operator L^^ : [H{n,V)]* —>■ 
H{fl,V) is well defined. Thanks to the compact embedding in lemma[l\i), we can consider the inverse operator L^^ 
restricted to L'^{fl) as a compact operator from L'^{Q,) into itself. 

Next, we will compare the eigenvalues defined on fl with the eigenvalues defined in the whole R^ . We will denote 
by H, D{L), cr„, 0„,- • • the space, the domain of the operator, the eigenvalues and the eigenfunctions and so on 
whenever the operator L is defined on the whole M^. More precisely, let us define the Hilbert space 



H := {u(^ iji(R^) : / V u^ < +oo 



and consider the operator L defined by 



Lu := -Ai 



Vu, 



for u e D{L) 



where now 



D{L) := Hn{u& L^{R^) : -Au + Vue ^^(M^)} 



(23) 



(24) 



(25) 



^), with compact embedding whenever N > 2 
admits a unique principal eigenvalue in M.^ 



Observe that the above lemmas [I] [2] still hold, in particular H ^^ L^{ 

and V satisfies hypothesis ^. The elliptic operator L as defined in 

noted by cti . This principal eigenvalue is the bottom of the spectrum of L in the function space H, and it admits an 

associated positive principal eigenfunction. 

In the following lemma we collect the monotonicity properties of the eigenvalues with respect to the domain, see |38l . 
Theorem 2.3] for a proof. As a consequence, we can compare the eigenvalues defined on fi C M^ with the eigenvalues 
defined in the whole M^. 
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Lemma 3. Let cri(ri) < cr2(^) < ■•• be the eigenvalues of L, with corresponding eigenf unctions (j)i{^), <f>2{^) ■ ■ ■ 
defined on H{fl). For a subdomain Q* d fl with boundary regular enough, let ai{fl*) < cr2(^*) < • • • be the eigenvalues 
of L, with corresponding eigenfuctions 4>i{^*)^ 4>2{^*) ■ ■ ■ defined on H(fl*), then 

a„(r!*) >(7„(r!). (26) 

In particular, let cti < o'2 < ■ ■ • be the eigenvalues of L, with corresponding eigenvectors (^„ defined on H, then 

an{n) > an- (27) 

The following lemma states the Maximum Principle for unbounded domains, see j4lL Theorem IX. 27]. 

Lemma 4. The maximum principle for the Dirichlet problem. Let il C M^ he an open nonempty set, possibly 
unbounded, with boundary of class C^, and assume V satisfy hypothesis 1^. Let f G L'^{fl) and u G H{i}, V) be such 
that m) holds. 
Then 

minjinf u, inf f}<u< maxjsup u, sup /} (28) 

on n Q^ Q 

where sup = sup ess and inf — inf ess. 
In particular, if 

u>0 on dVt, and / > in U, (29) 

then 

u>Oinn, and ||u||L~(f7) < niax{||M||i^(ao), ||/||l«=(o)}- (30) 

C. Main result 

We now consider Eq. ([6]) as a bifurcation problem, considering A as a real parameter, we look for pairs (A, u\) eM.xH 
such that u\ is a positive solution of ([5]). 

As in Sec. IV Bl let cti stand for the first eigenvalue of the eigenvalue problem 

--A + V{x)j(t>i:^ai<Pi, x G K^, (/)i(a;) ^ 0, as |a;| ^ oo, (31) 

and given an open regular enough domain fl C K.^, let ci = cri(ri) stand for the first eigenvalue of the Dirichlet 
eigenvalue problem 

-iA + F(a;)j 01 := 0-1(^1, x e n, (j)i{x) = 0, on dn, (32) 

where the first eigenfunction 0i = 0i (fi) . 

Let to the set where g vanishes, defined by Eq. ([5]). Technically we will require that this coefficient is bounded away 
from zero asymptotically, i. e. 

9{x) > 30 > 0, for |a;| > 1, (33) 

(note that in our numerical simulations we have observed the phenomenon described in the paper in less strict 
situations). In general, ui is not a connected set, let cti(u;) be the minimum 

aiiuj) := min{ai(wj) : 1 < j < J}, (34) 

where ujj : I < j < J are defined in (112). 

Theorem 1. The problem (0) has a positive solution (X,u\) whenever 

CTi < A < 0-1 (cj), (35) 

moreover 

\Wx\\l'^{rn)-^Q, as a i CTi, (36) 

II "A II L°- (K") ^ oo, as At CTi(w). (37) 
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Proof. The Crandall-Rabinowitz's bifurcation theorem, [42 implies that ((Ti,0) is a bifurcation point in M x if. 
Moreover, the bifurcation is supercritical, i.e. on the right of the eigenvalue. 
And, the monotonicity with respect to the parameter implies the uniqueness. 

A = (Ti(R^,y + guP-i) < (Ji{uj,V + guP-^) = a:i{uj,V). (38) 

D 

VI. CONCLUSIONS AND DISCUSSION 

We have studied the ground state of nonlinear Schrodingcr equations with confining potentials and repulsive inter- 
actions vanishing on nonempty sets lu and shown that there is a nonlinear localization phenomenon of the solution on 
this set. We have presented analytical formulae and numerical simulations showing it and provided a theorem which 
rigorously explains some of the observed features. 

The results have interesting implications to the physics of Bose-Einstein condensates with repulsive spatially varying 
interactions because in that situation the atom density will localize dramatically near the regions where the interactions 
are close to zero when the nonlinear interactions are relevant enough. This behavior does not depend on the spatial 
dimensionality and can be used to design very effectively one, two or three-dimensional spatial distributions with large 
values of the atom density by acting on the control field, e.g. using micromagnets to change locally the magnetic field 
or localized laser beams. 

The behavior described in this paper cannot be achieved with ordinary potentials since due to tunneling of the 
wavefunction in linear potentials it is not posinble to achieve such a high degree of localization. However, the nonlinear 
interactions here are essential in limiting the tunneling to the regions with larger values of the scattering length which 
is independent on the number of particles. 

Our results could be useful in different applications of Bose-Einstein condensates such as for atom lithography, 
atom beam guiding or other applications where a precise control of the positioning of large values of the density of 
an atomic cloud is required. 

These phenomena may appear in multicomponent condensates that offer a much wider range of possibilities for 
their ground states depending on the type of components and their relative interactions [47| . It would be interesting to 
study the situation in which the variations of the scattering length make the components to be miscible or inmiscible 
in different spatial regions and what is the geometry of the resulting domains. 

Finally, it is interesting that the phenomena reported here could also be observed in other systems ruled by nonlinear 
Schrodingcr equations similar to Eq. (|4]), for instance in nonlinear optical systems. Although ordinary optical materials 
have small nonlinearities were these phenomena could not be easily observed except for unrealistically large laser 
intensities, there are media with a giant nonlinear response [45l |46| as it happens when a probe laser beam propagates 
in a medium with transparency induced electromagnetically by a second coupling field. Therefore, managing the 
parameters it may be posible to find optical version of the phenomena presented here. 
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